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ABSTRACT 

The Reuven Ramaty High Energy Solar Spectroscopy Imager (RHESSI) has 
obtained the first high- resolution measurements of nuclear de-excitation lines pro- 
duced by energetic ions accelerated in a solar flare, a GOES X4.8 event occurring 
on 23 July, 2002 at a heliocentric angle of ~ 73°. Lines of neon, magnesium, 
silicon, iron, carbon, and oxygen were resolved for the first time. They exhibit 
Doppler redshifts of 0.1-0.8% and broadening of 0.1-2.1% (FWHM), generally 
decreasing with mass. The measured redshifts are larger than expected for a 
model of an interacting ion distribution isotropic in the downward hemisphere 
in a radial magnetic field. Possible interpretations of the large redshifts include 
1) an inclination of the loop magnetic field to the solar surface so that the ion 
distribution is oriented more directly away from the observer, and 2) extreme 
beaming of the ions downward along a magnetic field normal to the solar surface. 
Bulk downward motion of the plasma in which the accelerated ions interact can 
be ruled out. 

Subject headings: Sumflares — Sun:X-rays, gamma rays — line:profiles — gamma 
rays:observations 

1. Introduction 

We report the first high-energy-resolution measurements of nuclear de-excitation lines 
in a solar flare. These lines were detected by the Reuven Ramaty High Energy Solar Spec- 
troscopic Imager (RHESSI) during the GOES X4.8-class solar flare of 23 July, 2002, which 
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occurred at coordinates S13E72 on the solar disk (73° heliocentric angle). The characteristics 
of this flare are reviewed in depth in Lin et al. (2003). 

Gamma-ray lines from inelastic interactions of accelerated ions with ambient nuclei in 
solar flares (Lingenfelter & Ramaty 1967) were first observed by Chupp et al. (1973). The de- 
excitation lines can be Doppler shifted due to the nuclear recoils from the ion interaction and 
the emission of the gamma ray. The line profiles thus reveal the angular distribution of the 
interacting ions (Ramaty and Crannell 1976). Accelerated protons and a-particles produce 
narrow lines while accelerated heavy ions produce broad lines which merge together to form 
an underlying nuclear continuum. Share et al. (2002) analyzed an ensemble of gamma-ray 
line flares from the Solar Maximum Mission Gamma- Ray Spectrometer (SMM/GRS). They 
found that narrow lines exhibit ~ 1% redshift for flares at small heliocentric angles, but are 
not appreciably shifted near the limb of the Sun. Their measurements are consistent with a 
distribution of interacting particles which is isotropic in the downward hemisphere and zero 
in the upward hemisphere. 

Although the SMM/GRS data yielded useful measurements of some nuclear-line red- 
shifts and put constraints on the line widths, the energy resolution of that instrument was 
moderate, averaging around 4% FWHM over the nuclear-line range. RHESSI (Lin et al. 
2002) was designed to provide high-resolution measurements of the shape and redshift of 
de-excitation lines using cryogenically cooled germanium detectors to achieve an energy res- 
olution averaging around 0.2% FWHM from 1-6 MeV (Smith et al. 2002). 

2. Analysis and Results 

The RHESSI spectrometer consists of nine coaxial germanium detectors, each divided 
into a ~1 cm upper segment (which stops most hard x-rays) and a ~7 cm rear segment. 
Here we use only data from the rear segments, which have most of the efficiency for stopping 
gamma-rays above 300 keV (Smith et al. 2002). Detector #2, which is operated in an un- 
segmented mode and has degraded energy resolution, was not used. The spectra in Figure 1 
were accumulated from 00:27:20 to 00:43:20 UT on 23 July, 2002, an interval that includes 
most of the high-energy emission from this flare (Lin et al. 2003). 

The flare intensity is roughly half the background level above ~1 MeV. The instrumen- 
tal gamma-ray background is caused by the interactions of cosmic rays and protons trapped 
in the radiation belts with materials in the spacecraft and the atmosphere. Every 15 or- 
bits (about one day), the geomagnetic coordinates of the spacecraft approximately repeat 
themselves, and therefore so do the cosmic and trapped particle fluxes. Thus we model the 
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background during the flare with the average of two spectra, taken 15 orbits before and after 
the flare, each of which is the same length as the flare observation itself. Above 3 MeV, how- 
ever, we replaced this background spectrum with a different one using a full day's integrated 
data. In this high-energy band the background is simple in form and nearly free of lines, 
so that we can can get an accurate background estimate with small statistical fluctuations, 
regardless of the exact orbital parameters, simply by rescaling the full-day integration to 
match the intensity of the ± 15-orbit spectrum. 

The full spectrum (shown in Lin et al. (2003)) from 250 keV to 8.5 MeV was fitted 
with six bright de-excitation lines (Table 1) plus a broken power-law representing electron 
bremsstrahlung (the index hardens from 2.77 to 2.23 at a break energy of 617 keV). It 
was also necessary to represent two other components of the spectrum: the very broad 
lines from interactions of accelerated heavy ions with ambient hydrogen, and the faint lines 
from other proton- and a-induced lines that cannot be measured individually. Within the 
statistical limitations of the data, the sum of these components is well-represented by three 
very broad Gaussians centered at 1845, 4358 and 6575 keV. Bright narrow lines were included 
from positron annihilation (Share et al. 2003a) and neutron capture (Murphy et al. 2003) 
at 511 keV and 2223 keV, respectively. Finally, weak lines at 1263 keV and 6918 keV of 
intermediate width were also included, representing expected features which are compounds 
of more than one line. 

The fitting was done with the Spectral Executive (SPEX) package (Schwartz 1996), 
taking into account the full diagonal and off-diagonal response of the instrument (Smith et 
al. 2002). The spectra for the six bright de-excitation lines are shown in Figure 1 and the fit 
parameters are given in Table 1. The features at ~840 keV in the upper left panel of Figure 1 
and ~1807 keV in the lower left panel may be the 844 keV and 1809 keV solar de-excitation 
lines of aluminum (Kozlovsky, Murphy & Ramaty 2002), but if so are at a surprisingly 
high flux. The effect of the instrumental resolution (about 3.0 keV at 847 keV, 4.1 keV at 
1779 keV, and ~ 10 keV at 6129 keV) has been removed from the width measurements, 
but, except for the 56 Fe line, it is negligible compared to the solar line widths. The best-fit 
Gaussians and their zero-redshift versions are also shown in Figure 1. Zero redshift is visibly 
inconsistent with the data for most lines. 

The redshifts and widths from Table 1 are shown in Figure 2. Heavier nuclei will 
recoil less from a collision with a proton or a-particle, and will thus show less redshift and 
broadening. This trend is clearly visible. The expected relation is not exactly linear, since 
the average energy of the accelerated particles interacting with each nucleus depends on the 
the cross section for excitation. The neon line, which appears to have a redshift below the 
trend in Figure 2, has the lowest threshold for excitation. 
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Fig. 1.— RHESSI background-subtracted count spectra from 00:27:20 UT to 00:43:20 UT 
on 23 July, 2002. Each panel is labeled with the element primarily responsible for the line 
shown. The carbon and oxygen lines also show the secondary peak from escape of a 511 keV 
positron-annihilation photon, which also contains information on the line shape. The heavy 
curve shown in each panel is the Gaussian fit from Table 1 plus the underlying bremsstrahlung 
continuum and broad lines (see text), convolved with the instrument response. The lighter 
line is the same fit forced to zero redshift for comparison. The error bars are ±lcr from 
Poisson statistics. 
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Fig. 2. — Redshift (top) and Doppler broadening (bottom) from Table 1 plotted as a function 
of the reciprocal of the mass of the nucleus. Data points marked as squares are from this 
flare. Data points marked with dashes are the average of five flares observed by SMM/GKS 
near 74° heliocentric angle. 
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3. Discussion 

The sum of the fluences in the narrow lines of Table 1 is 137 ph cm -2 . Share & 
Murphy (1995) calculated the total fluence in narrow lines for 19 X-class flares observed 
by SMM/GRS. Estimating the contribution of fainter lines not in Table 1 by comparison 
with the flares studied by Share & Murphy (1995) gives a total line fluence of about 180 ph 
cm -2 for this flare, similar to the brighter events from SMM/GRS, which include flares from 
GOES class X2.8 to X15.0. 

Table 1 and Figure 2 also show SMM/GRS results from Share et al. (2002) for the average 
of five flares close to a heliocentric angle of 74°. The most striking difference is the lack of 
redshift in the lines from carbon and oxygen, in significant disagreement with our result for 
the 23 July 2002 flare. The null redshift and high broadening in the SMM/GRS neon line 
are in marginal disagreement with RHESSI as well. The advantage of high energy resolution 
is clear: the RHESSI results are more precise using a single flare than the SMM/GRS result 
for the average of five flares, even though RHESSI is a smaller instrument. With more 
flares observed at high resolution, we will be able to tell whether this disagreement between 
RHESSI and SMM/GRS is due to intrinsic variations between flares or to the limitations of 
SMM's moderate-resolution spectroscopy. 

We use the modeling code of Murphy, Kozlovsky & Ramaty (1988), which allows an 
arbitrary angular distribution of interacting particles about an axis at an arbitrary angle 
to the observer (9). For particles distributed with respect to a radial field line, 9 is equal 
to the heliocentric angle, 73° for this flare. The simulation generates separate gamma-ray 
spectra for a-particles and protons, and includes the effects of nuclear recoil both from the 
interaction and the emission of the gamma ray. Thus, even a purely downward beam of 
incident protons gives broadened lines due to the range of recoil angles. 

We simulated two angular distributions for the interacting ions: 1) isotropic within the 
forward hemisphere and zero in the hemisphere toward the observer, and 2) beamed directly 
forwards. Since we will explore the possibility that the magnetic field is not perpendicular 
to the solar surface, we will use "forward" to mean along the magnetic field away from the 
observer and "downward" to mean perpendicular to the solar surface. 

Distributions with no forward bias, such as a fan beam or a full isotropic distribution, 
give no redshift and thus disagree strongly with the redshift values in Table 1. For each of the 
distributions, the viewing angle 9 was run at 0°, 30°, 45°, 52°, 60°, and 73°. Interpolations 
along cos((9) gave line shapes for intermediate angles. The proton and a-particle energy 
spectra were modeled as power laws with index -3.75, and the a/proton ratio was fixed at 
0.5, both consistent with results from SMM/GRS (Share & Murphy 1995, 1998). Figure 3 
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Fig. 3. — Model 20 Ne line shape for a power-law index of -3.75, a viewing angle of 30°, and 
a forward-isotropic distribution. Dashed line: line shape from interacting protons. Dash- 
dotted line: line shape from interacting a-particles. Solid curve: total line shape for an 
a/proton ratio of 0.5. Vertical line: the rest energy of the de-excitation line, 1634 keV. 
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shows the 20 Ne line shape predicted by the forward-isotropic model for a viewing angle of 
30°. 

For both models and each viewing angle, we fitted the model lineshapes (e.g. Figure 3) 
to the data, allowing the line ffuences and the underlying broad lines and continuum to vary 
to minimize x 2 ■ Figure 4 shows the relative x 2 f° r the two angular distributions and all 
viewing angles. The forward isotropic model is clearly a much better fit, but it fits best 
with a viewing angle of 30-40°. If the field lines of the loop were perpendicular to the solar 
surface at their base, where the ions interact, the viewing angle would be the heliocentric 
angle of the flare, 73°. This is inconsistent with the data at greater than 99% confidence. 
RHESSI's ability to image gamma rays allows us to be certain that the line emission comes 
from near the flare, i.e. near 73° (Hurford et al. 2003). 

We have run further simulations with power-law indices from -1.75 to -4.75 and a/proton 
ratios from 0.0 to 2.5. In all cases the best-fit viewing angle for the forward isotropic case is 
close to 30-40° and inconsistent with 73°. Note that even the forward beam, which produces 
the maximum possible redshift of any distribution, still fits best with a viewing angle smaller 
than the heliocentric angle. For spectral index -2.75 and harder, the forward beam can fit 
as well as the forward isotropic distribution. 

A redshift could be caused either by recoil from an anisotropic distribution of interacting 
ions or by bulk motion away from the observer of the medium in which they interact. We 
believe the latter is not significant here for two reasons. First, in the case of bulk motion 
the percentage redshift is a constant, and not inversely proportional to the mass of the 
excited nucleus (Figure 2). Second, although downward flows in the chromosphere have 
been observed in impulsive flare events via red wings in Ha line profiles, the velocities are 
only 20-50 km/s (e.g. Wiilser, Canfield & Zarro 1992), implying redshifts of no more than 
0.007-0.017%, a trivial fraction of the total gamma-ray shifts. 

Two explanations remain for the large redshifts: either the field lines in the loop legs 
are inclined toward us instead of being perpendicular to the solar surface, or else the angular 
distribution of the interacting particles is closer to a pure forward beam than a forward- 
isotropic distribution. Forward-isotropic or isotropic distributions of interacting a-particles 
are strongly favored over a downward beam in a study of the lines at 0.429 and 0.478 MeV 
from a-a reactions in this flare (Share et al. 2003b) and in flares observed with SMM/GRS 
(Share & Murphy 1997). Also, the low average redshifts of the 20 Ne, 12 C and 16 lines in all 
the flares seen by SMM/GRS (Share et al. 2002), including those used in Table 1, favored 
a downward isotropic distribution (or one with a slight additional upward component) to a 
downward beam. 
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Fig. 4. — Relative values of x 2 f° r model fits to the 23 July flare data at various viewing 
angles. Solid curve: forward-isotropic distribution. Dashed curve: forward beam. Dotted 
lines: confidence levels for the hypothesis that the best fit is preferred to a given model. The 
left edge of the plot corresponds to a loop viewed side-on, the right edge to a loop viewed 
down along the legs, and the vertical line corresponds to a loop perpendicular to the local 
solar surface. 
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If, on the other hand, magnetic loops containing accelerated ions are inclined from the 
perpendicular for this flare, it may help explain why the SMM/GKS widths are larger than 
ours in Table 1. Share et al. (2002) combined five flares to get their result near 74°, and if 
these had a wider range of 9 than of heliocentric angle, differing redshifts could combine to 
make a broader line. Solar active regions often do contain inclined loops, and the magnetic 
topology of regions such as this one ("P^fS") that produce the largest flares also tend to be 
the most complex. 

A more detailed analysis of nuclear de-excitation line shapes can constrain not only 
the angular distribution of the particles but also the a/proton ratio (as can be seen from 
the different component shapes in Figure 3) and even the spectral index. The full power 
of the spectral analysis of gamma-ray lines, however, will be realized when the constraints 
imposed by the line shapes are combined with information from the gamma-ray line fluences 
and from other observations such as gamma-ray and x-ray imaging (Hurford et al. 2003), 
magnetograms and imaging at other wavelengths, and the decay profile of the neutron- 
capture line (Murphy et al. 2003). 

The work at the University of California, Berkeley and at NASA's Goddard Space Flight 
Center was supported by NASA contract NAS5-98033, and that at the Naval Research 
Laboratory by NASA DPR W19746. We thank Hugh Hudson, Gordon Hurford, Brian 
Dennis, George Fisher, and Benzion Kozlovsky for useful discussions. 
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Table 1. Best fit Gaussian parameters for prompt nuclear lines 



Isotope 


Rest Energy 
(keV) 


Fit Energy 
(keV) 


% Redshift 


FWHM 
(keV) 


% FWHM 


Fluence 
(ph cm -2 ) 


SMM/GRS 
% Redshift a 


SMM/GRS 
% FWHM a 


56 Fe 


847 


846.091°;™ 


u - 1J --0.07 


1 2+ 2 ' 9 


14+ - 34 
u - i4 -0.13 


7 5 +3.4 
' -°-2.3 






24 Mg 


1369 


1363.6^20 


n 40+ - 17 

U - 4U -0.14 


2l.0±|;0 


1 54+O.59 


OQ q + 7.2 

2 °-:>-6.6 






20 Ne 


1634 


1628. 8lH 


S2+ ' 10 




1 07+ ' 26 
J- 1 " -0.22 


21.4±|;| 


0.0 ±0.2 


2.9 ±1.0 


28 Si 


1779 


1776.812;? 


12+ ' 11 


16 7+ 4 ' 5 
1D - '-5.4 


n 04+0.25 
u - a4 -0.30 


17.111" 






12 C 


4438 


44031 1° 


fl 7Q+0-23 
u - ' a -0.22 


QO+ 42 
3 -29 


2 06+ ' 95 


28.6li 3 6 1 


-0.07±0.14 b 


3.18 ±0.74 


16Q 


6129 


6094lJ| 


U.OO_ 2 g 


122±|f 


1 QQ+ 111 
1 - aa -0.83 


34.21J 2 ; 8 


-0.26 ± 0.13 b 





a Average of five flares near heliocentric angle 74° 
b Blueshift 



